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ABSTRACT 

This paper, the second of a series on radiation-regulated accretion onto black holes (BHs) 
from galactic scales, focuses on the effects of radiation pressure and angular momentum of 
the accreting gas. We simulate accretion onto intermediate-mass black holes, but we derive 
general scaling relationships that are solutions of the Bondi problem with radiation feedback 
valid for any mass of the BH Mbh- Thermal pressure of the ionized sphere around the BH 
regulates the accretion rate producing periodic and short-lived luminosity bursts. We find that 
for ambient gas densities exceeding oc M^^, the period of the oscillations decreases 

rapidly and the duty cycle increases from 6%, in agreement with observations of the fraction of 
active galactic nuclei at z ~ 3, to 50%. The mean accretion rate becomes Eddington limited for 
nn^oo > '^h'Ix) — ^Hoo-^oo^ where Too,a is the gas temperature in units of 10^ K. In the sub- 
Eddington regime, the mean accretion rate onto BHs is about 1%)T^\ of the Bondi rate, thus 
is proportional to the thermal pressure of the ambient medium. The period of the oscillations 
coincides with depletion timescale of the gas inside the ionized bubble surrounding the BH. Gas 
depletion is dominated by a pressure gradient pushing the gas outward if nH,oo < ''^h oo 
accretion onto the BH otherwise. Generally, for nH,oo < '^h oo angular momentum does not 
affect significantly the accretion rate and period of the oscillations. 

Subject headings: accretion, accretion disks - black hole physics - dark ages, reionization, 
first stars - hydrodynamics - methods: numerical - radiative transfer 



1. INTRODUCTION 



Gravitational accretion onto point sources can be described analytically ( Bondi & Hoyle|1944 Bondi 



1952 1 assuming spherical symmetry. The Eddington-limited Bondi formula is often prescribed in cosmo- 



logical simulations to estimate the gas accretion rate onto black holes (BHs) from large scales ( [Volonteri & 



Rees|2005t[Di Matteo, Colberg, Springel, Hemquist, & Sijacki|2008t|Pelupessy, Pi Matteo, & Ciardi|2007t 



Greif, Johnson, Klessen, & Bromm|2008l [Alvarez, Wise, & Abel|2009l |Kim, Wise, Alvarez, & Abel|201l] . 



However, even in sub-Eddington regime, a fraction of the gravitational potential energy of accreted gas is 
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converted into mechanical or radiative feedback ( [Shapiro|1973 1, reducing the accretion rate. The radiation 



emitted by BHs creates feedback loops that regulate the gas accretion and luminosity of the BHs. Several 
published works have investigated physical processes that may dominate the feedback such as X-ray preheat- 
ing, gas coohng, photo-heating, and radiation pressures ( jOstriker, Weaver, Yahil, & McCray||1976[ |Covyie, 



Ostriker, & Stark||1978; Bisnovatyi-Kogan & Blinnikov||1980 


Krolik & London|1983[ Vitello||1984 Wan- 


del, Yahil, & Milgrom|1984;|Milosavljevic, Bromm, Couch, & Oh|2009a 


Ostriker, Choi, Ciotti, Novak, & 



Proga 2010 ; Novak, Ostriker, & Ciotti|20TT] |. In general, radiative feedback reduces the accretion luminosity 



of the accreting BH ( [Ostriker, Weaver, Yahil, & McCray|197^ |Begelman|1985[ [Ricotti, Ostriker, & Mack 



2008] ) ■ There have been extensive publications on self -regulation of supermassive BH growth at the centers 



of eUiptical galaxies (^Sazonov, Ostriker, Ciotti, & Sunyaev|2005| 


Ciotti & OstrLker|2007 Ciotti, Ostriker, & 


Proga||2009; Lusso & Ciotti||2011 1 or axisymmetric outflows in active galactic nuclei (AGNs; 


Proga||2007 


Proga, Ostriker, & Kurosawa 2008 ; Kurosawa, Proga, & Nagamine 


2009 


Kurosawa & Proga 


2009a|bl. 



Recently, several works have paid closer attention to radiation-regulated accretion onto intermediate-mass 
black holes (IMBHs; p"ilosavljevic, Bromm, Couch, & Oh|20"09at [Milosavljevic, Couch, & Bromm|2009b 



Park & Ricotti|2011t |Lil20TT] ). 



Cosmological simulations show that massive BHs may have formed in metal-free minihalos as Pop- 



ulation III star remnants in the early universe (Abel, Anninos, Norman, & Zhang 1998 Bromm, Coppi, 
& Larson||1999t [Abel, Bryan, & Norman||2000t [Madau & Rees||200l| [Schneider, Ferrara, Natarajan, "& 
Omukai|2002[ |0h & Haiman[[2002| | or from direct collapse of primordial gas ( [Carr, Bond, & Amett|1984 



Haehnelt, Natarajan, & Rees[[I998[ [Fryer, Woosley, & Heger[|2001t [Begelman, Volonteri, & Rees[[2006 



Volonteri, Lodato, & Natarajan|2008t[Omukai, Schneider, & Haiman|2008t [Regan & Haehnelt[2009t [Mayer," 
Kazantzidis, Escala, & Callegari|2010} Johnson, Khochfar, Greif, & Durier'2011). Estimating the accretion 
luminosity of IMBHs (for a review, see I Miller & Colbert|2004 van der Marel|2004| ) is important to under- 
stand their cosmological importance at high z and in the local universe ( Mack, Ostriker, & Ricotti[[2007t 
Ricotti|2()09 l. Since the luminosity of IMBHs is directly related to their accretion rate, these studies are also 



relevant for better understanding the mass growth of primordial massive BHs in the early universe (Madau 
& Rees[[200T] [Volonteri, Haardt, & Mad"au[[2003t [Yoo & Miralda-Escude[[2004t [Volonteri & Rees|[2005l 
Johnson & Bromm[[2007t [Pelupessy, Pi Matteo, & Ciardi[[2007l [Alvarez, Wise, & Abel[[2009| ) or provide 
clues about the origin and impact on the ionization history of the universe of ultraluminous X-ray sources 
(ULXs; [Krolik, McKee, & Tarter[[T98l] [Krolik & Kailman[[T984t [KroUk[[2004t [Ricotti & Osttiker|[2004l 
Ricotti, Ostriker, & Gnedin|2005t[Ricotti|2007[[Strohmayer & Mushotzky|2009[ ). 



In Park & Ricotti (201 1 1, hereafter Paper I, we explored the radiation-regulated accretion onto IMBHs 
assuming spherical symmetry and zero angular momentum of the accreting gas. One of the main objectives 
of the study was to derive an analytical description of the effect of radiation feedback on the Bondi accretion 
rate. We accomplished this goal by simulating accretion onto IMBHs with idealized initial conditions and 
simple physics, but for a large parameter space of the initial conditions (varying BH radiative efficiency, 
the BH mass, the density and temperature of ambient gas, and the spectrum of radiation). We found that 
the IMBH is quiescent most of the time with short intense periodic bursts of accretion (with duty cycle 
of 6%). The qualitative description of the cycle is as follows. Gas accumulates in a dense shell ahead of 



-3- 



the I-front nearly halting accretion onto the IMBH. Meanwhile, the ionized gas inside the hot bubble is 
pushed outward toward the dense shell by pressure gradients, eventually de -pressurizing the hot bubble, 
producing the collapse of the shell and a burst of accretion. The scaling relationships for the burst period, 
mean and peak accretion rates can be understood analytically but are quite sensitive to the details of the 
thermal structure inside the Stromgren sphere. Thus, we expect that gas metalUcity may be an important 
parameter in the problem that we have not yet explored. 

In this paper, the second of the series, we relax most of the simplifying assumptions in Paper I and dis- 
cuss the effects of helium heating/cooling, radiation pressure, and gas angular momentum on the accretion 
rate. It has been noted that not only electron scattering but also radiation pressure on H I may be important 
( Milosavljevic et al.|2009a i. We explore how the radiation pressure regulates the gas accretion by transfer- 



ring momentum to the inflowing gas, and whether these physical processes become important compared to 
the pressure gradients inside the Stromgren sphere. 

As in Paper I of this series, here we present the results of our simulations as dimensionless accretion 
rates Arad = M /Mb, where Mb = Tre^/^pooC^ M^j^c^^ is the Bondi accretion rate for an isothermal gas 
(with polytropic index 7 = 1), that is a function of BH mass Mbh, density p^o, and sound speed of c^^oo of 
neighboring gas. The Eddington luminosity is L-^dd = 47rGMbh?TT^pCcrr^^ , which is proportional to the BH 
mass only. We define the Eddington accretion rate Mgad = -^Edd/c^ which is a factor of 10 smaller than 
the definition most often used Mgdd = L^dd/O-lc^, in which rj = 0.1 is assumed, and the dimensionless 
Eddington rate as Agdd = ^~^-^Edd/-^s- We introduced three quantities to describe the accretion: mean 
accretion rate (Arad), accretion rate at peaks Arad,max» and period between accretion bursts Tcyde- In Paper I 
we found 

(\ad) ^ C(nH,oo) [Tin(a, Z)/4 x 10^ K] , (1) 
with C ~ 3% for nn.oo > 10^ cm'^ and C ~ 3%(nH,oo/10^ cm-3)i/2 f^j. ^^^^ ^ '^q5 ^.^-3 ^j^^ 7] = 0.1 
and Mbh = 100 Mq. Here, Tin{a, Z) is the time-averaged temperature at the accretion radius within the 
Stromgren sphere, that is 4 x lO'' K for our fiducial spectrum (spectral slope a = 1.5) and a gas composed of 
hydrogen only (metallicity Z = 0). However, Tin is sensitive to changes of a and the gas composition (see 

Paper I and Section 3 in this paper). In the sub-Eddington regime, the period Tcydo oc {Rg) where {Rs) is 

1/2 

the time-averaged Stromgren radius, and the duty cycle /duty = Ton/Tcycie = (Arad)/Arad,max ~ 6% T^' 4, 
where Ton is the duration of bursts. 

This paper is organized as follows. In Section 2, we briefly explain the numerical methods. In Section 3, 
we present our simulation results including the effect of helium cooling/heating, radiation pressure, and gas 
angular momentum. Summary and discussion are presented in Section 4. 



2. NUMERICAL SIMULATIONS 



We run a suite of hydrodynamic simulations to interpret how radiative feedback regulates accretion 
onto BHs. We use a modified parallel non-relativistic hydrodynamics code, ZEUS-MP (Stone & Norman 
1992[ Hayes et al.|2006[ ) plus a radiative transfer algorithm ( Ricotti et al.|2001 1 to simulate photo-ionization 
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and photo-heating by UV and X-ray ionizing photons emitted near the BHs. See Paper I for a detailed 
description. 

In this study, we include the effects of helium heating/cooing in addition to hydrogen. Therefore, we 
simulate photo-ionization, photo-heating and coohng for six species H I, H II, He I, He II, He III, and e~. 

We also calculate the radiation pressures both on and H I to interpret the effect of momentum 
transfer to the inflowing gas by the ionizing photons. The magnitude of acceleration at a given radius due to 
radiation pressure depends on the luminosity, the ionization fraction of hydrogen and helium, and the cross 
section of the species to photon-ionization. The specific flux F,^ oc e^^/r^ at a given radius(r), assuming 
a power-law spectrum with a spectral index a, depends on the optical depth Tj^, and the cross section ajy. 
Thus, the accelerations due to momentum transfer to H l and can be written as 

arad,H I = / <7H \MFydu, (2) 



arad,e- = / CFrFudv, (3) 

TUpC J 

where xh i and Xg- are H I and e~ fractions, respectively, (Tt is the Thomson cross section, and nip is the 
proton mass. The radial component of the acceleration at a given radius is updated as a = agrav + a^ad. 
where a^ad = arad.H i + aj-ad.e- ■ 



In Section 3.3 we study the effect of non-zero angular momentum of gas which leads to a time delay 
between the accretion rate at the sonic radius and the luminosity output, due to the formation of an accretion 
disk. In order to estimate realistic values of the time delay we assume that the gas conserves angular 
momentum and settles into an accretion disk of radius i?disk- We then assume an alpha model for the thin 
disk to estimate the timescale for the gas to lose angular momentum and fall into the BH. 

Numerically, it is convenient to express the time delay in units of the free-fall timescale calculated at 
the simulation's inner boundary (typically i?min ^ 10^^ pc). The free-fall timescale we have defined can be 
very large compared to calculated at the radius of the accretion disk near the BH (at tens of gravitational 
radii R^ch = 2GMbh/c^). Approximately, the gas is accreted at the viscous timescale tvisc, that compared 
to % is 

tvisc(i?disk)/tff(i?disk) ~ a'^M^ ~ ""^c^^.^kG'Mbhiidisk ~ 0-5 a"nc/cs,disk)^^disk' (4) 



where a is the dimensionless parameter for a thin disk ( [Shakura & Sunyaev 1973 1, Cs,disk is the sound speed 



of the gas in the disk, and we define T^disk = ^disk/^Sch- The dependence of the free-fall time on radius 
is ifj oc R^-^, while the viscous timescales as tvisc oc -R^^tfr oc R^'^ assuming constant sound speed due to 
effective cooling (note that since we are considering a gas of zero or very low metallicity, the gas in the disk 
will not easily cool to temperature below 10^ K if the gas is atomic). Thus, the infall time at the disk radius 

-Rdisk is 

tvisc(^disk) 0.5 t4in 1/2 0.3 / Tdisk \ ^1/2 
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To estimate the parameters in Equation (jsjl we have defined t^min = (GMbh/^min)^''^ — 260 km s ^. 
Assuming a ~ 0.01-0.1, Tdisk ~ 10"^ K, and IZd^sk < lO^-lO"^, we find time delays of < 300 free-fall 
times at i?min. that is the parameter space we explore in Section [33] 

In our code, the accretion rates calculated at the inner boundary of the simulations are stored in an 
assigned array about 1000 steps for each Stored accretion rates with a given time delay are then read 
from the array and used to estimate the luminosity at the current moment. 



3. RESULTS 

In this section, we show the results of simulations discussing the effects of helium heating/cooling, 
radiation pressures, and angular momentum on the BH accretion rate. All the simulations in this paper 
include helium heating/cooling but the gas is metal free. The simulations remain qualitatively the same as 
in Paper I where we did not include helium; the only noticeable difference with respect to Paper I is that 
the accretion rate at peak luminosity shows multiple minor peaks instead of a well-defined single peak. 
This is to be expected, as a larger opacity produces a stronger feedback with respect to the hydrogen only 
case, leading to multiple shocks in the gas. This complicate structure - i.e., a burst consisting of several 



sub-bursts - is commonly found (e.g., Ciotti & Ostriker 2007 1. In addition, the average accretion rate (Arad 



decreases from ~ 3% to ~ 1%, but this can be understood by the increase of the mean temperature inside 
the Stromgren sphere to from Tin ~ 4 x 10^ K to ~ 6 x 10"^ K. The top panel of Figure[T]shows the accretion 
rate as a function of t/ = 0.01-0.1 with Mbh = 100 Mq, nn,oo = 10^ cm'^, and Too = 10^ K. Large 
symbols show (Arad) while small symbols show Arad,max- For the given set of parameters, the luminosity 
remains in the sub-Eddington regime, thus the effects of radiation pressures are minor. The bottom panel of 
Figure [l] shows the dependence of Tcydc on rf^/'^, the same as found in Paper I. However, Tcycic for r] = 0.1 
is now ~ 2200 years which is ~ 60% of the value found in Paper I for the given set of parameters. This is 
also well understood (see Equation (22) in Paper I) as our model predicts Tcyde oc (Arad)^^^- 



3.1. Effect of Radiation Pressures 

In Paper I, we have focused on exploring the parameter space in which the mean accretion rate is 
dominated by thermal feedback, i.e., radiation pressure can be neglected. We found (Arad) ~ 1% for 
'^H.oo = 10^ cm~^, assuming Mbh = 100 Mq, T^o = 10^ K, a = 1.5, and including helium cool- 
ing/heating. However, not surprisingly, including the effect of radiation pressure produces a reduction of 
the accretion rate when the BH luminosity approaches the Eddington limit. Figure [2] shows (Arad) as a 
function of gas density for a 100 Mq BH, comparing simulations that do not include radiation pressure 
(open triangles) to ones including pressure on H I only (open squares), on e~ only (open pentagons), and 
the total effect of radiation pressure (solid circles). Compton radiation pressure reduces the accretion rate 
below (Arad) ~ 1% for nH,oo ^ 10^ cm~^ while the radiation pressure on H I appears always negligible 
with respect to Compton scattering. Both (Arad) and Arad.max change from a constant fraction of the Bondi 



-6- 



accretion rate to the Eddington rate AEdd> shown by the dashed Une for Mbh = 100 Mq and radiative 
efficiency r] = 0.1. 

Figure [3] shows the dimensionless accretion rates (Aj-ad) and Arad,max as a function of the BH mass 
from Mbh = 10^ to 10^ Mq, keeping the other parameters constant: r] = 0.1, nH,oo = 10^ cm^^, and 
Too = 10*^ K. The simulations include radiation pressures on H I and e^, and show that the transition to 
Eddington-limited accretion happens for Afbh — 5000 Mq. 



3.1.1. Transition from Bondi-like to Eddington-limited Accretion 

So far the simulation results have shown that Compton scattering on electrons is the dominant radiation 
pressure effect, thus the Eddington-limit applies. Figure|4]summarizes the results of a large set of simulations 
that include radiation pressure. The top three panels in Figure ^a) shows (A^ad) as a function of gas 
density for Mbh = 10^, 10'^, and 10^ Mq, respectively. For each BH mass, corresponding Eddington 
limits are shown by the dashed lines. The panels show the mean accretion rate (Aj-ad) (large triangles) and 
Arad,max (small triangles) transitioning from being a constant fraction of the Bondi rate at low densities to 
being Eddington-limited at higher densities. The period of the accretion Tcydo, in the bottom panels, also 
shows different dependencies in Bondi-like and Eddington-limited regimes. We will come back to this in 
Section 1X21 

Figure Qb) shows the mean accretion luminosity in units of Lsdd for Mbh = 10^, 10^, and 10^ Mq 
as a function of gas density. The dotted lines show 1% of the Bondi accretion rate for each BH mass. Thus, 
from Figure |4] approximately we have 

(M) = min(l%r^5^MB, //-^MEdd), (6) 

where Too,4 = Too/(10^ K), valid for density nn.oo ^ 10^ cm^'^, and a = 1.5. 

It is thus apparent that IMBHs can grow at a rate near the Eddington limit if the gas density of the 
environment is larger than the critical density 

nH,^~4xl0cm ^^^j j (-) . (7) 



3.1.2. Why is Continuum Radiation Pressure Negligible ? 

As shown in Figure 2]|4 the simulations show that radiation pressure on H I does not play an im- 
portant role when the accretion rate is sub-Eddington. In this section, we focus on understanding why 
this is. Figure [5] shows the evolution of relative magnitude of acceleration due to radiation pressures 
normalized by the gravitational acceleration at a given radius. Each panel refers to a different density 
"-H.oo = 10^, 10^, 10^, and 10® cm^'^. Within the Stromgren sphere, the relative effect of Compton radi- 
ation pressure remains constant as a function of the radius since the electron fraction Xg- is close to unity 
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and the gas is nearly transparent to ionizing radiation. Outside of the Stromgren sphere, the rapid decrease 
of the electron fraction reduces the effect of Compton scattering. Radiation pressure on H I (thick lines in 
Figure |5]) increases as a function of radius and has its peak value just inside the Stromgren sphere. This is 
due to the increase of the H I fraction as a function of radius. Outside the Stromgren sphere the relative 
effect of H l radiation pressure drops quickly because the ionizing luminosity decreases rapidly due to the 
increase of the H i opacity. 

Continuum radiation pressure on H I is comparable to Compton electron scattering only in a shell just 
inside the Stromgren sphere, where the H I abundance starts to increase rapidly as a function of radius 
and the ionizing radiation is not fully shielded by H I. With increasing gas density, the peak and mean 
luminosities increase, hence the relative effect of Compton pressure on average increases and eventually 
becomes comparable to the effect of gravity (i.e., Eddington limit). Whereas Figure[5]shows that the relative 
effect of continuum radiation pressure does not increase much with increasing gas density. In addition, the 
range of variation of radiation pressures during a period of oscillation decreases with increasing density. 
In other words, at low densities radiation pressures display several magnitudes of variation which are not 
seen in the high-density regime. As a result, at low densities (nu,oo ^ 10^ cm~'^) radiation pressure 
is significant only near the peaks of luminosity and generally is negligible compared to gravity; whereas 
at high densities (nH,oo ^ 10^ cm~'^) Compton scattering dominates throughout a period of oscillation 
reducing the accretion rate to Eddington-limited values. Only at intermediate densities nH,oo — 10^ cm~^, 
the magnitude of H I radiation pressure just behind the Stromgren radius becomes comparable to that by 
Compton scattering. 

The weak dependence of the H l radiation pressure on density and its magnitude with respect to the 
Compton pressure can be understood analytically. The key point is that the H I radiation pressure is propor- 
tional to the value of the neutral fraction xu i just behind the Stromgren radius Rg and, assuming ionization 
equilibrium, it is easy to show that xh i{Rs) oc n^^^^. It follows that the pressure on H I is relatively 
insensitive to variations of nn: 

PSd « Soxu i{Rs) exp [-t{Rs)] oc n^/^ (8) 

where 5*0 oc nn is the ionizing luminosity, and exp [—t{Rs)] = const. The derivation of xn i{Rs) is as 
follows. At Rs = SQ^^n-J^^^a^^^^^ the photoionization rate is T{Rs) = Soa^jAirR'^ oc n^^. Assuming 
photoionization equilibrium xh i{Rs)^{Rs) = nuaji, we demonstrate that 

fry X rinaR -2/3 
XH i{Rs) = N oc ' . (9) 



3.2. Two Self-regulated Modes of Accretion: Collapsing I-front versus Quasi-steady I-front 

One of the most interesting aspects of the radiation-regulated accretion onto BHs is the qualitative 
change of the period and duty cycle of the luminosity bursts observed in the high-density regime. As argued 
in Paper I and confirmed by further simulations in this work, the physical reason for this transition is a 
change of the dominant mechanism depleting the gas inside the Stromgren sphere between two consecutive 
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bursts. In the low-density regime, gas is pushed outward toward the ionization front by a pressure gradient 
(hereafter, mode-I accretion). At higher-densities gas accretion onto the BH becomes the dominant gas 
depletion mechanism (hereafter, mode-II accretion). Incidentally, as discussed below, simulations show 
that radiation pressure becomes important near the transition from mode-I to mode-II, at least for most 
of the initial conditions we have simulated. In Paper I, we have observed mode-II accretion only for our 
highest density simulation (for nn.oo = 10^ cm~'^ and M^h = 100 Mq). In this paper, to better understand 
this regime, we have extended the parameter space to higher densities and higher BH masses. Figure [6] 
shows snapshots of the density (top halves in each panel) and ionization fraction (bottom halves in each 
panel) for two-dimensional simulations including radiation pressure, for nH,oo = 10^ cm^"^ (top panels) 
and nH,oo = 10"^ cm^^ (bottom panels). The snapshots are taken for each simulation at the moment of 
a burst of the accretion rate (left panels), in-between two bursts (middle panels), and just before a burst 
(right panels). For ambient density nH,oo = 10^ cm~^, the Stromgren sphere collapses onto the BH which 
leads to a strong luminosity burst. On the contrary, the size of Stromgren sphere does not change much 
during the oscillation period for ambient density nn.oo = 10^ cm~^. In this latter case, the oscillation of the 
accretion luminosity is driven by density and pressure waves originating at the I-front, while in the former 
case, the collapse of the I-front onto the BH leads to a much more intense accretion burst. In Figure [TJleft), 
we compare the accretion rate onto the BH as a function of time for nH,oo = 10^ cm^'^ (top panel) and 
'^H.oo = 10^ cm^^ (bottom panel). For nH,oo = 10^ cm^'^, the collapse of I-front leads to strong burst 
of gas accretion, with Arad,max about x20(Arad)- Hence, the duty cycle /^^^^ = (Arad)/Arad,max is about 
6%. The pressure gradient inside Stromgren sphere supports the gas shell accumulating at the I-front from 
collapsing until the accretion rate drops 4-5 orders of magnitude compared to the accretion during the burst. 
However, the Stromgren radius remains remarkably constant before its collapse due to the decline of gas 
density inside the H II region. In contrast, in the nH,oo = 10"^ cm^'^ simulation the accretion rate peaks at a 
few times (Arad) before decreasing by about 2 orders of magnitude. The duty cycle approaches f^^^^ ~ 50% 
for this mode of accretion. As shown in Figures [3] and [4j simulations that do not include radiation pressure 
also show a rapid decrease of the period Tcydc and Arad,max with increasing gas density and BH mass, but 
the mean accretion rate (Arad) does not. Thus, the reduced value of Arad,max/ (Arad) = l//duty explains the 
longer duty cycle observed for mode-II accretion. A more detailed illustration of the qualitative difference 
between mode-I and mode-II accretion is shown in Figure [TJright). The figure shows the time evolution 
of the gas density profile (top panels), the temperature profile (middle panels) and the hydrogen ionization 
fraction (bottom panels) for the reH,oo = 

cm and njj co — 10 cm simulations. Small variations 
of the density, temperature, and ionization fraction profiles are observed for nn.oo = 10^ cm^'^, while clear 
collapses of I-front are observed in the evolution of the profiles for tiu qq = 10 cm Note that this 
quasi-stationary profile is not produced by the effects of radiation pressures. The same effect is found for 
'^H.oo = 10'' cm~^ without including radiation pressure effects. 

Interestingly, for our fiducial case simulations (Mbh = 100 Mq, = 10^ K, = 0.1, and a = 1.5), 
the critical density at which the mean accretion rate becomes Eddington-limited nearly coincides with the 
critical density for transition from mode-I to mode-II accretion. This explains why the mean accretion rate 
and the peak accretion rate become Eddington-limited at nearly the same density. Indeed, if while increasing 
nH,oo> the duty cycle remained at about 6% as in mode-I accretion, the mean accretion rate would not be 
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able to approach the Eddington limit, even though the peak accretion can be mildly super-Eddington. We 
will show below that the transition to mode-II accretion depends on the free parameters in the problem and 
may take place at much lower densities than the critical density for Eddington-limited accretion. 



derstand why there exists a clear transition to the Eddington-limited regime with increasing density or BH 
mass. For mode-I accretion, radiation pressure may become comparable to the gravity near the Stromgren 
radius, but this effect dominates only for a short time, during the peaks of luminosity. The peak accretion 
can indeed become moderately super-Eddington for a short time, also because of the broken spherical sym- 
metry of the collapsing shell due to Rayleigh-Taylor instability of the accreting gas. However, for mode-II 
accretion, the geometry of accretion from large scales is quasi-spherical and radiation pressure effects are 
significant during the most of the duration of oscillations, hence the accretion rate is Eddington limited. 

Figure [8] shows the relationship between the period of accretion bursts and the average size of the 
Stromgren sphere produced by the accreting BH. When the gas depletion inside the H II region is dominated 
by the outward flow of gas toward the I-front, t^^^.^^ shows a linear relation with {Rg) (solid line). This 
linear relation is almost identical to the results in Paper I, where helium cooling/heating was not included. 
By increasing the ambient gas density, eventually the gas depletion becomes dominated by accretion onto the 
BH. In this latter case, assuming that the dimensionless accretion rate (Arad) is constant (a valid assumption 
in the sub-Eddington regime), Tcydc scales as (Rs)^ (dotted line). However, the simulation results at high 
ambient gas density shown in Figure [s] are not well fitted by Tcyde oc (Rs)^, and indeed seem to follow a 
linear relationship Tcydc oc {Rs), similar to the low-density one but with an offset. This can be explained 
because at high densities the accretion rate becomes Eddington-limited soon after the transition to mode-II 
accretion for which Tcyde {Rs}"^- It follows that the assumption (Aiad) ~ const becomes invalid and 



instead T^yde = ^in « p{Rs)^ / MEdd- In this regime, since the Stromgren radius is (Rs)'^ oc r/MEdd/p^, we 
getpoc A4£(i?,)-i-5and 



for different values of Mbh- 

Thus, the small offset in Tcydc observed in Figure [8] when the density is increased, can be understood 
because ^, at which the transition from mode-I to mode-II accretion takes place, is nearly equal to nj^'^^, 
the critical density at which the mean accretion rate becomes Eddington limited. But, in general, the ratio 
of these critical densities may depend on all the free parameters of the model. 

Our analytical model of feedback-regulated feeding of the BH, can help understand the dependence 
of the critical density on all the parameter space, not fully covered by the simulations. We found that 
the cycle period Tcyde is the shortest time between the gas depletion timescales tin = Mn u/M, where 
-^H II ~ Pm{Rs)^ is the mass inside the H II region, and tout ~ 3(i?s)/cs,in (see Paper I). Thus, by 
definition, when the density approaches the critical density we have tin — ^out, but this condition also 
implies that the mean Stromgren radius approaches the effective accretion radius: 



The quasi-stationary I-front observed for the nn,. 



10^ 



cm 



^ simulation is also important to un- 




(Rs) 



cr 



10 X rb,eff- 



(11) 
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Equation (jllj), is derived setting M = (M) = 47r/3inCs,in?"b eff = (-^rad)-^^B in the relationship for tin- Since 
in our model we have AnTln - Poo^'oo, it follows that rb,eff ~ (71n/Too)^/^(Arad)^''^»^b, and {Rg)'''' ~ 
2T^^4Ti~'^J^rb, where = GM/cl^ is the Bondi radius, Too,4 = Too/lO^ K, and Tin,* = Ti^/Q x 10'' K 
is the mean temperature at the accretion radius inside the H II region (normalized to the value found for 
a = 1.5). Thus, {RsY^ and period of the bursts are 

{R,r ^ (0.01 pc)Mbh,2Ti-y', (12) 
r-,!, ^ (1000yr)Mbh,2Ti-'/', (13) 

with Afbh,2 = Mbh/100 Mq. 

Applying naively the analytical expression for the Stromgren radius produced by a source of luminosity 
L = rjc'^{M) in a gas of density poo gives {Rs) oc M'^^^n-^^^T^^rj^/^. However, using the simulation data, 
we find that the mean radius of the Stromgren sphere in the sub-Eddington regime is nearly independent of 
Too, and if ~ (i?s)'^'^, is also independent of rj: 

(0.015 pc)Mja(^)-"^(^)-"\ 

where E = Lq/Sqis the mean energy of ionizing photons, and we have assumed hydrogen recombination 
coefficient aR = (4 x lO^"' cm^ /s)T^^^]/'^ . In addition, we find that, (Rs) oc r]^^^ as expected for {Rs) S> 
{Rsy^. The deviation from the naive expectation is not surprising, as the BH luminosity and the density 
inside the Stromgren sphere are not constant with time. Indeed, although both the maximum and mean 
luminosities of the BH are oc r], the simulations show that the luminosity at the minimum of the cycle, L™™, 
is nearly independent of t]. Typically ^^^''^ <^ L, but when ?2h,oo 

approaches the critical value L™™ ~ L. 
Similarly, assuming an effective mean density (pmPoo)^^^ in the Stromgren radius expression would explain 
the temperature dependence in Equation ([14]). 



Finally, setting (Rs) = {Rs)'^'^ we derive the critical density 



ng,oo ~ (5 X lO*^ cm-3)M-,;,r£: ' . (15) 

The critical density as well as other scaling relationships in our model depends on E = L / So 
and Tin, but for a gas of zero metallicity (including helium), these quantities are determined only by the 
spectrum of the radiation. Assuming a power law spectrum with index a is easy to show that 



a/ {a — 1) if a > 1, 



E = 13.6 eV < 



ln(/ii^max/13.6 eV) ifa = l (16) 

a/(l-a)(/ifmax/13.6eV)" ifa<l. 



We have estimated /iz^max = 0.2 keV as the frequency at which the mean free path of the photons equals 
{Rs). The points in Figured show Tin as a function of E for simulations with a = 0.5, 1, 1.5, 2, 2.5 taken 
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from Figure 9 in Paper I. The line shows the fit to the points: 



E -^"^ 



For our fiducial model, for which £' ~ 41 eV, the value of the critical density is very close to given 
in Equation Q: 

H,oo ^ / 

From an inspection of Equation ( 1 8 1 is evident that the only cases in which ng ^ can be larger than 
are assuming the largest realistic values of unity for Too,4 and and assuming a spectrum of radiation 
from the BH softer than a = 1.5 that would reduce E with respect to the fiducial value. Vice versa a 
hard spectrum, low radiative efficiency and accretion from a gas colder than 10^ K would decrease the ratio 
'^H oo/'^hIx) below unity, making mode-II accretion sub-Eddington for a wider range of densities. For these 
cases the period of the cycle could become very short with increasing density as Tcyde {Rs)'^ oc n^^^- 



3.3. Effect of Non-zero Angular Momentum of Gas 

As discussed in Paper I, the introduction of small angular momentum in the flow, which is realistic in 
most astrophysical problems, can modify the time-dependent behavior of accretion rate presented in this se- 
ries of papers. Angular momentum of gas leads to the formation of an accretion disk near the Schwartzshild 
radius of a BH. This disk is not resolved in our simulations. Thus, the accreted gas may experience a time 
delay before it is converted to radiation. Here, we test how the introduction of time delay would affect the 
feedback loops of accretion. 

As mentioned in Section 2, it is important to estimate physically motivated time delays. Here, we 
explore the time delay of 1-300 times tfr(i?min) which is large enough with an assumption of a-disk model. 
On the other hand, no matter how long is the time delay, what really matters is how the time delay compares 
to the oscillation period, which depends mainly on the gas density for a fixed mass of BH. We investigate this 
issue in the low (nn.oo = 10^ cm^^) and high density (nn.oo = 10^ cm^'^) regimes where the oscillation 
pattern and the periods are different. At low densities a time delay of a few hundred free-fall times is much 
smaller compared to the oscillation period, whereas at high densities the maximum time delay that we have 



tested is comparable to the oscillation period. In the left panel of Figure 10 which shows the result for 
'^H.oo = 10^ cm~^, Tcycie docs not increase at all as a function of time delay since the introduced time delay 
is much smaller than the original oscillation period. In the right panel of Figure 10 for nH,oo = 10"^ cm~^. 



the maximum time delay that we introduce is comparable to the original oscillation period, and we see that 
Tcycie increase approximately by the amount of time delay. In both cases, we still observe oscillations. Thus, 
only in the case of accretion from a high-density gas which produces shorter oscillation period, and for an 
accretion disk with i?disk ~ ^min> the time delay may have an important effect on the accretion rate. 



Indeed, the accretion disk may not only introduce a time delay but also smooth out the accretion rate 
on a timescale of the order of the viscous timescale. In this case for cases in which the disk is large (i?disk ~ 
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i?niin) and Tcycio is short (i.e., for mode-II accretion), the disk may further smooth out or completely erase 
the periodic low-amplitude oscillations in the accretion rate from large scales. 



4. Summary and Discussion: Scaling Relationships 

We have presented a systematic study on how the classic Bondi problem of spherical accretion onto 
a compact object is modified by the effects of radiation feedback. We solve radiative transfer equations in 
the radial direction for the hydrogen and helium ionizing radiation emitted by the BH. Gas is optically thin 
inside Stromgren radius while it becomes optically thick for gas outside the ionized gas. In this paper, the 
second of a series, we have focused on the effects that radiation pressure and angular momentum have on 
the gas supply and accretion rate onto the BH. The simulations focused on accretion onto IMBHs but the 
analytical scaling relationships we have derived are rather general, and although the initial conditions are 
somewhat idealized, should describe reality more accurately than the classical Bondi formulae. 

Here, we summarize the main results and scaling relationships we found in the first two papers of this 
series for non-moving BHs accreting from a uniform medium. In our models we have assumed that the BH 
accretes from a uniform density and temperature reservoir, significantly larger than the Bondi radius and 
{Rs)- This assumption is well motivated for accretion onto stellar and IMBH, but for SMBH there could be 
supply of gas to the BH from stars within (Rs) (stellar winds) or other astrophysical object (merger-driven 
accretion, etc). The scaling relationships can be applied to problems involving a wide range of masses of 
the accretor, from stellar mass objects (e.g.. Wheeler & Johnson||2011 1 to supermassive BHs. One caveat 



is that we are neglecting the effects of self-gravity of the gas (see Li||201 1 1 and the gravitational potential 



due to the dark matter halo of the host galaxy, which may play an important role for the case of accretion 
onto supermassive BHs. Indeed, a simple calculation shows that at the I-front gravity due to the mass of 
the gas inside the Stromgren sphere exceeds the BH's gravity if Mbh ^ 10^ MQ/(r/_iToo,4). Our model 
predicts scaling relationships for the period, duty cycle, peak and mean accretion onto the BH, as well 
as relevant critical densities and size of the Stromgren sphere around the BH. In the following summary 



of the scaling relationships, we express Tin in the equations in terms of E given by Equation ( 17 1 that is 



valid for a gas of low-metallicity. E is related to the spectral index a by Equation ( 16 1. For higher values 
of the gas metallicity, the coefficients in the equations can be different due to changes in the relationship 
between and the spectrum of the radiation. A caveat is that our simulations have explored a large but 
limited parameter space for the masses of the BHs, temperature and density of the ambient gas, etc. So, the 
proposed scaling relationships, although they are based on a physically motivated model we inferred from 
the simulations, should be used with caution for sets of parameters that are significantly different from the 
range confirmed by simulations. 



The main qualitative result of our study is that radiation feedback produces periodic oscillations of the 
accretion rate from large scales onto the BH, and thus periodic short-lived bursts of the BH luminosity. We 
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found two modes of self -regulated accretion, determined by 



5 X 10^ cm 

Mbh,2 



E 



41 eV 



-9/16 



(19) 



If '^H.oo < oo (mode-I), the accretion luminosity of the BH has regular bursts with period T^fy^j^ during 
which the BH increases its brightness by about 5 orders of magnitude but only for a short fraction of the cycle 

T 1 /2 

period: the duty cycle is 7^^^^ = Ton/Tcycic ~ 6% Tj^ ^. During the quiescent phase in the accretion cycle 
the gas accumulates in a dense shell in front of the H II region rather than accreting directly onto the BH. As 
the luminosity decreases after the burst, the density inside the H II region also decreases because is pushed 
outward by a pressure gradient, thus maintaining the I-front radius nearly constant. Eventually the density 
and pressure inside the H II region cannot sustain the weight of the dense shell that collapses producing a 
burst of accretion. The cycle repeats regularly. If nn.oo > ^hoo (mode-II) the cycle is qualitatively different: 
the duty cycle is about /^Jj^y > 50% and the peak accretion rate is only a few times the mean. There is no 
collapse phase of the dense shell and the H II region remains roughly stationary while the accretion rate 
oscillates. The physical motivation for mode-II accretion is that the timescale for the depletion of the gas 
inside the H II region becomes dominated by accretion onto the BH. Only for mode-II accretion the BH 
growth rate can approach the Eddington limit, given that the density exceeds the critical density 



n 



Edd 
H,oo 



4 X 10^ cm 



rr-1 -1 



(20) 



For nearly all realistic cases > 



n 



H,oo- 



For Mbh = 100 Mq, at densities 10^ cm"'^ < nH,oo < "hIx) mean accretion rate onto the 
BH, in units of the Bondi rate is (Arad) ~ '^%T^\{E/4:1 eV)~^, independent of all the other parameters. 

For nH,oo < 10^ cm^'^ instead, (Arad) ~ 1%("-H,oo/I0 



cm 



) 1/2^^5^ (£:/41 eY)-\ depends weakly 
on the gas density. One caveat is that in Paper I the dependence on the free parameters of the transition 
density nH,oo = 10^ cm~^ has been only partially explored. As shown in FigureHof the present paper, the 



simulation results are consistent with a transition density inversely proportional to the BH mass. Hence, if 
5 A/f-i ^,^-3 ^ ^„ < ^hIxd mean accretion rate is proportional to the thermal pressure nn.ooT'oo 



lO^iVV;^ cm 



of the ambient gas: 



(M)«(4xl0i8gs-i)M2 



105 



cm 



^oo,4 



E 



41 eV 



(21) 



If nH,oo > ng'i^ then (M) = L^ddivc"^)'^ ■ The duty cycle is 



duty 



./Lv^50% ifnH,oo>nf 



/duty ~ ^"-"^ "H.oo / "-H.oo' 

and the maximum accretion luminosity which depends on the duty cycle thus is, 



(22) 



IT 



mm 



Edd 



l,^7?_lMbh,2 



105 



cm" 



T ■ ^ 



(23) 
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where LEdd = 1-3 x 10'*'^Mbh,2 erg s ^, and 



\^"«0.06 ifnH,oo>nfo^ 
The cycle of the oscillations also falls into two regimes: 

^ ^ f-cVcie-(0.1Myr)M,^/5,^f(^)-V3(_^ ifnH,oo<n-^ 
^'^'^ V^cle-(lGy^)r/-l(aF^)-H4I^v)-'/^ if-H,oo>ng,^. 

The astrophysical apphcations of this model are innumerable and are beyond the aim of this paper to 
discuss them in detail. However, one of the most obvious results is that the luminosity of an accreting BH 
should be smaller than the value inferred applying the Bondi formula. Not only because the mean accretion 
rate is always < 1% of the Bondi rate, but also because if nH,oo < "-h oc ^^^^ about 5 orders 

of magnitude lower than the Bondi rate inferred from the ambient medium temperature and density. Thus, 
this simple arguments could have interesting consequences to interpret the observed quiescence of SMBH 

1/2 

in ellipticals and Sgr A*. Also, the duty cycle of ~ 6% ^ we found for mode-I accretion is interestingly 



close to the fraction of galaxies with AGNs ~ 3% found deep field surveys (e.g., Steidel et"aL]|2002 Luo 
. For SMBHs of about 10^ Mq, ng^ ~ 500 cm-^. 



etal. 



2011 



In this paper we also found that IMBH can grow at near the Eddington limit if nH,oo > ™ax(n^'^j^ , 
This has potentially important consequences on the ability of seed IMBH from Population III stars to grow 
by accretion into SMBH during the first gigayears of the universe age. This possibility seemed precluded if 
the duty cycle of the burst was 6% as found in previous works. 

Finally, although the nature of ULXs is unknown, there are indications that they may host an IMBH 
(e.g., jStrohmayer & Mushotzky||2009 1. An IMBH accreting from an interstellar medium (ISM) with high 



pressure such as dense molecular cloud (nH.oo^oo ~ 10^-10"^ cm^'^ K) would be L™^^ ~ lO'^'^-lO'^^ erg 
for Mbh = 1000 Mq, that is comparable to the luminosity of ULXs. However, this assumes that the IMBH 
is at rest with respect to the ISM. We will focus on gas accretion onto moving BHs with radiation feedback 
in the third paper of this series (K. Park & M. Ricotti 2012, in preparation). Clearly, more work is needed 
to address each of the aforementioned topics in detail, but the basic ground work presented in the present 
paper may allow the re-visitation of a few longstanding problems still unsolved in astrophysics. 



The authors thank the referee, Luca Ciotti, Richard Mushotzky, Chris Reynolds, Eve Ostriker, and 
Edward Shaya for constructive comments and feedback. The simulations presented in this paper were carried 
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Technology at the University of Maryland ("deepthought"). This research was supported by NASA grants 
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Fig. 1 . — Top: accretion rate in units of the Bondi rate as a function of radiative efficiency rj for simulations 
with Mbh = 100 Mq, f^H,oo — 10*^ cm and — lO'^ K. Large symbols indicate mean accretion 
rate ((Arad) ~ while small symbols show accretion rate at peaks (Arad,max ~ 20%). Bottom: period 
between bursts Tcydc as a function rj. The dotted line shows Tcydc oc 77^/^. In both panels, triangles represent 
simulations neglecting the effect of radiation pressures, while circles show simulations including radiation 
pressures. Radiation pressures introduce a minor difference in both the accretion rate and period of the 
bursts for this parameter set. 
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Fig. 2. — Comparison of relative importance of radiation pressures in regulating mean accretion rate (top 
panel) and the period of bursts (bottom panel) as a function of the ambient gas density nH,oo- Symbols 
are explained in the legend of the figure. Top: long dashed line shows the Eddington limit for a 100 Mq 
BH with r] = 0.1. When the accretion rate is sub-Eddington ( nH,oo < 10® cm^^ ) radiation pressures 
both on electrons and H I do not play an important role and the thermal structure of the Stromgren sphere 
regulates the accretion. Radiation pressure is important in reducing the accretion rate at nu,oo = 10"^ cm^^ 
where the accretion rate approaches the Eddington rate. The transition of accretion rate from (Arad) ~ 1% 
to the Eddington-limited regime happens at ~ 4 x 10® cm^'^ for a 100 Mq BH with r/ = 0.1 and 

Too = 10^ K. Bottom: radiation pressures do not produce significant differences in Tcyde- Transition of 
Tcycie from mode-I (dotted line) to mode-II (short dashed line) happens at the critical density ng ^ ~ '^h*!^ 



(see Section 3.2 1. The result shows a good agreement with the work of Milosavljevic et al. (2009b I 
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Fig. 3. — Same as Figure [TJ but showing (Arad)> -^rad,max» and Tcycic as a function of Mbh with r] = 0.1, 
'^H.oo = 10^ cm~^, and Too = 10^ K. A similar pattern which we observe as a function of density is also 
seen as a function of Mbh- With increasing Mbh, the transition from (Arad) ~ 1% to the Eddington-limited 
regime and the transition of Tcyde from mode-I to mode-II happen at Mbh ~ 4 x 10^ Mq. 
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Fig. 4. — Left: same as Figure |2j but showing (Arad)> Arad,max> and Tcyde as a function of gas density 
for Mbh = 10^ Mq (left panel), 10^ Mq (middle panel), and 10^ Mq (right panel). Long dashed hnes 
in each panel show the Eddington accretion rate for rj = 0.1 and the given BH mass. With increasing gas 
density, the accretion rate eventually becomes Eddington limited, but the transition to the Eddington-limited 
regime occurs at densities ~ 4 x 10^ cm^M^^ decreasing linearly with increasing BH mass. Right: 
accretion luminosities normalized by Eddington luminosities for the same simulations in the left figure. 
Symbols (circles: 10^ Mq, triangles: 10^ Mq, squares: 10^ Mq) show the simulations including radiation 
pressures for each BH mass. With increasing gas density, the accretion rate becomes regulated primarily by 
Compton radiation pressure. 
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log (radius(pc)) 

Fig. 5. — Radial profiles of the gas acceleration due to radiation pressures on H I and e~ normalized 
to the gravitational acceleration of simulations for tih.oo = 10^10^10^and 10^ cm with BH mass 
Mbh = 100 Mq, r] = 0.1, and Too = 10^ K. Thick lines refer to radiation pressure on H I, while thin 
lines show Compton scattering radiation pressure. Different line types show the profiles at different time 
during the oscillation cycle (e.g., solid lines at the accretion bursts and dot-dashed just before the bursts). 
Radiation pressure on H I peaks just inside the Stromgren sphere with weak dependence on density. While 
Compton radiation pressure inside the Stromgren sphere increases on average as a function of density. At 
'^H.oo = 10' cm^'^ the peak values of H I radiation pressure and Compton radiation pressure become 
comparable and about 10% of gravity. However, at higher densities (> 10^ cm~^) radiation pressure on 
electrons becomes dominant everywhere inside the Stromgren sphere. 
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Fig. 6. — Illustration of the two different modes of oscillations found at different ambient densities. The 
panels show the time evolution of gas density and ionization fraction in 2D simulation for a BH of mass 
Mbh = 100 Mq, gas density nH,oo = 10^ cm~'^ (top panels), and nH,oo = 10*^ cm~^ (bottom panels). 
In each panel, top halves show the density and the bottom halves show the ionization fraction, xh n = 
ii/^H- When the density is nH,oo ^ oo — 5 x 10^ cm^, the collapse of ionization front onto the BH 
leads to a burst of accretion luminosity (mode-I). For densities 71h,oo > '^H.oo the size of Stromgren sphere 
does not change much with time (mode-II). Note the different oscillation modes of the accretion rate and 
luminosity are driven by the collapse of dense shell (mode-I) and by a density wave (mode-II). 
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Fig. 7. — Left: accretion rates as a function of time for nn.oo = 10^ cm~^ and nH,oo = 10^ cm~'^ assuming 
r] = 0.1, Mbh = 100 Mq, and Too = 10^ K. Different modes of oscillations occur at different density 
regimes. Mode-I oscillation at nH,oo = 10^ cm~^ shows about 5 orders of magnitude range between peak 
and the minimum accretion rate, while mode-II oscillation at nH,oo = 10^ cm~'^ shows only 2 orders of 
magnitude range. Right: evolution of radial profiles for density (top panel), temperature (middle panel), 
and neutral/ionization fractions (bottom panel) of the same simulations in the left figure. Note the change 
of physical properties inside Stromgren sphere during a period of mode-I oscillation (nH,oo = 10^ cm^'^), 
while mild changes are observed for mode-II oscillation (nH,oo = 10^ cm~^). 
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Fig. 8. — Relationship between the period of the accretion bursts, Tcycie, and the time-averaged size of 
Stromgren radius {Rg). Tcyde shows a Unear relation with (Rg) when the gas depletion inside the Stromgren 
sphere is dominated by a pressure gradient inside the H II region that push the gas toward the I-front. 
Instead, Tcydc {Rs}"^ (dotted lines for each ilfbh) when the gas depletion is dominated by accretion onto 
the BH. With increasing density of the ambient gas, for each Mbh, the transition to mode-11 oscillation 
and the transition to Eddington-limited regime happen at similar densities nH,oo = "-§00^ ^h'Ix.- 
Eddington-limited regime Tcyde becomes proportional to (i?^)^/^ for each BH mass (dashed lines). 
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Fig. 9. — Temperature Tin at the effective inner Bondi radius, located inside the H ll region produced by the 
accreting BH, as a function of the mean energy of ionizing photons E of the spectrum of radiation emitted 
near the BH by the accretion disk. We have assumed a gas of nearly zero-metallicity and a power-law 
spectrum Fj, oc 



10 



t /t 

delay' cycle.O 

0.0001 0.001 0.01 0.1 0.01 0.1 



10 



1 — I I I 1 1 1 ll 1 — I I I 1 1 1 ll 1 — I I I 1 1 III 1 — I I I 1 1 1 1 



I I I 1 1 1 ll 1 — I I I 1 1 III 1 — I I I 1 1 1 ll 1 — I I I 1 1 III 



n„„ = 10^ 



n„„ = lO' cm"'^ 



I 



t 



Q 2 I I I I I I I I I ll I I I I I I I l_ 

1 10 100 1000 1 10 100 1000 



^delay/^ff.l 



Rmin 



Fig. 10. — Dependence of the period between bursts Tcycie on the time delay between the accretion rate 
at -Rmin and the BH output luminosity. The time delay is produced by the presence of the accretion disk 
that is unresolved in our simulations. The two panels show Tcydc in units of rcycie,o = ^cycic (*deiay = 0) 
as a function of tdciay for nH,oo = 10^ cm~^ (left panel) and 10^ cm~^ (right panel). The bottom axis 
shows ideiay in units of the free-fall time at Rmin and the top axis in units of rcycie,o- The introduction 
of a time delay does not change Xcyde when the gas density is nu,oo = 10^ cm~^, while Tcycie increases 
by approximately the amount of time delay introduced for nu,oo = 10^ cm~^. In this density regime, the 
largest time delays introduced are comparable to the oscillation period rcycie,o- In both cases, the oscillatory 
behavior of the accretion luminosity does not disappear. 



